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Abstract 

This paper deals with a dynamical system that generalizes the Kepler- Coulomb system and 
the Hartmann system. It is shown that the Schrodinger equation for this generalized Kepler- 
Coulomb system can be separated in prolate spheroidal coordinates. The coefficients of the 
interbasis expansions between three bases (spherical, parabolic and spheroidal) are studied in 
detail. It is found that the coefficients for the expansion of the parabolic basis in terms of the 
spherical basis, and vice-versa, can be expressed through the Clebsch-Gordan coefficients for 
the group SU(2) analytically continued to real values of their arguments. The coefficients for 
the expansions of the spheroidal basis in terms of the spherical and parabolic bases are proved 
to satisfy three-term recursion relations. 
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1 . Introduction 



The purpose of the present paper is to further study the nonrelativistic quantum mechanical 
system corresponding to the three-dimensional axially symmetric potential 

^ ® / •> i 5 i 9 + ft / 9 i ') i 9 9 i 9 + ^ 9 I 9 ("0 

y/x z + y z + z z \/x l + y z + z z x z + y z x z + y z 
(in Cartesian coordinates), where a, (3 and 7 are constants such that a < and 7 > If 
/3 = 7 = 0, we have the ordinary (spherically symmetric) Kepler-Coulomb potential. In the 
case where (3 = and 7 7^ 0, the potential (1) reduces to the Hartmann potential that has been 
used for describing axially symmetric systems like ring-shaped molecules [1-3] and investigated 
from different points of view in the last decade [4-20] . The generalized Kepler-Coulomb system 
corresponding to (3 7^ and 7 7^ has been worked out in Refs. [9, 14, 17-20]. In particular, 
the (quantum mechanical) discrete spectrum for the generalized Kepler- Coulomb system is well 
known [9, 14, 17], even for the so-called (q, p)-analogue of this system [17]. Furthermore, a path 
integral treatment of the potential (1) has been given in Refs. [14, 18]. Recently, the dynamical 
symmetry of the generalized Kepler- Coulomb system has been studied in Refs. [17, 19, 20], the 
classical motion of a particle moving in the potential (1) has been considered in Ref. [19], and 
the coefficients connecting the parabolic and spherical bases have been identified in Ref. [20] as 
Clebsch-Gordan coefficients of the pseudo-unitary group SU(1, 1). 

The potential V [see Eq. (1)] belongs to a set of three-dimensional potentials systematically 
investigated in the 1960s (in connection with the question of accidental degeneracy) by the late 
Professor Smorodinsky and some of his students and/or collaborators [21-23] and revived, in 
recent years, by Evans [24-26]. 

The generalized Kepler-Coulomb system turns out to be interesting from both a classical- 
mechanical and a quantum-mechanical point of view. Indeed, not all the (bounded) classical 
trajectories are periodic and the quantum energy spectrum exhibits accidental degeneracies 
with respect to the geometrical group 0(2). (The latter two points are probably connected.) 
The accidental degeneracy arises from the fact that the Schrodinger equation for the potential 
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V is separable in more than one system of coordinates [21-23]. In addition, the potential V 
generalizes the Kepler-Coulomb potential, that is of central importance in quantum chemistry, 
and, like the Hartmann potential, may have applications [possibly in a (q, p)-deformed form or 
in a supersymmetric form] in chemical physics for systems presenting a line of singularity. 

A further motivation for the present work is the following. It has been recognized for a long 
time that the Schrddinger equation and the Hamilton- Jacobi equation for the potential (1) can 
be separated in spherical polar and parabolic rotational coordinates [23]. However, to the best 
of our knowledge, it has never been questioned, even in the recent papers of Refs. [9, 14, 17- 
20], if these equations are also separable in an additional coordinate system as the system of 
prolated ellipsoidal and hyperboloidal coordinates. This is actually the case and the originality 
of our work is mainly concerned with this separability in prolate spheroidal coordinates. 

The system of spheroidal coordinates constitute a natural system for investigating many 
problems in mathematical physics (see Ref. ]27j and references cited therein). In quantum 
mechanics, the spheroidal coordinates play an important role because they are appropriate in 
describing the behavior of a charged particle in the field of two Coulomb centres. The distance 
R between the centres is a dimensional parameter characterizing the spheroidal coordinates. 
These coordinates are changed into spherical and parabolic coordinates as R — > and R — > oo, 
respectively, if the positions of one Coulomb centre and of the charged particle are fixed when 
taking the limits. In this sense, the spheroidal coordinates are more general than the spherical 
and parabolic ones. This explains the interest of spheroidal coordinates in various domains 
(astrophysics, plasma physics, theory of the chemical bond, etc.). Such an interest is well 
documented in quantum chemistry, for instance in the study of two-center systems like Hjj" 
and H2 (see Refs. [28-34] for a nonhexhaustive list of papers showing the interest of spheroidal 
coordinates in quantum chemistry). 

The paper is organized as follows. In Sections 2 and 4, we describe the spherical and 
parabolic bases, for the generalized Kepler- Coulomb system, in a way adapted to the introduc- 
tion of interbasis expansions. In Section 3, we prove an additional orthogonality property for 
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the spherical radial wavefunctions of given orbital quantum number I. In Section 5, by using 
the property of bi-orthogonality of the spherical basis, we calculate the coefficients of interbasis 
expansions between spherical and parabolic bases. In Section 6, we construct the spheroidal 
basis of the potential (1) in terms of the superposition of the spherical and parabolic bases. 
Finally, we prove in the appendix that the bi-orthogonality property (discussed in Section 3) of 
the spherical radial wavefunctions is a consequence of the accidental degeneracy of the discrete 
spectrum for the generalized Kepler-Coulomb system. 

2 . Spherical Basis 

By introducing the nonnegative constants q (i = 1, 2) such that 

P = c 2 - ci, 7 = ci + c 2 , (2) 
the potential V [see Eq. (1)] can be rewritten in the Poschl- Teller form 

in spherical coordinates (x = r sin # cosy?, y = r sin 9 sin f, z = rcos6). The Schrodinger 
equation 

= Em (4) 
for the potential (3) may be solved by searching for a wavefunction in the form 

*(r,e,<p) = R(r)Z(8,<p). (5) 

This leads to the two coupled differential equations 

l d ( . dz\ l d 2 z ( Cl c 2 \ 



1 d f 2 dR\ A 
r 2 dr \ dr J r 2 



;R + 2(e-^)R = 0, (7) 

where A is a separation constant in spherical coordinates. [We use a system of units for which 
the reduced mass ji and the Planck constant h satisfy ji = h/(2n) = 1.] 
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The solution of Eq. (6) is easily found to be 



Z{0,<p) = Z lm (e,i P ;S 1 ,5 2 ) = ^a5i,52)(l + cos^^(l-cose)^?y' |m|+il) (cose)e^, 

(8) 

where I G N, m G Z, 5i = \J m 2 + 4q — \m\ (for i = 1,2) and P^ 01 '^ denotes a Jacobi polynomial. 
Furthermore, the separation constant A is quantized as 



x.(, + *±*)(, + *±* + i). 

The normalization constant Ni m (8i,S 2 ) in (8) is given (up to a phase factor) by 

✓ . m— Itt 

-1 2 



(9) 



iV Im ((Ji,<J 2 ) 



21' 



A 



(2/ + 5i + 6 2 + - \m\)\T(l + \m\ + 5 X + 5 2 + 1) 



(10) 



2 «5i+fe+2 r (/ + 5l + !) r (/ + 52 + !) 

The angular wavefunctions Z; m [see Eq. (8)] shall be referred to as ring-shaped functions. These 
ring-shaped functions generalize the functions studied by Hartmann [1-3] in the case (3 = (i.e., 



Si = 5 2 ). Due to the connecting formula [35 



A + £) C*(s) = (2\) n p( X - hX ~^ (x) 



(11) 



between the Jacobi polynomial p( a >^ and the Gegenbauer polynomial C£, the case 5i = S 2 = 5 
yields 



Zi m (0,ip;5, 5) 



m — m 



-l)^^2 |m|+ T m \+S 



(21 + 25 H 




m )! 


4vr 2 r(/ + 


m 


+ 25+1) 



(sm0)H+<5d™^( C( ,^ )( >''^ 



'|_| m | ^<joi/;c , (12) 

a result already obtained in Ref. |36|]. [In Eq. (11), (x) n stands for a Pochhammer symbol.] The 
case 5 = (i.e., (3 = 7 = 0) can be treated by using the connecting formula 



r H + - (i-^)-i-c 



m| + i 
-H ' 



(13) 



- \ 2, 

between the Gegenbauer polynomial and the associated Legendre function P ; ' m ' [ jJJJ . In fact 
for 5 = 0, Eq. (12) can be reduced to 



1 

,y;0,o) = (-1)- 



21 + 1 (/ - |m|)! 
V 4tt (/ + |m|)! M 



P m| (cos^e 



irrup 



(14) 



an expression that coincides with the usual (surface) spherical harmonic Yi m (9,ip) (e.g., see 



[37|). The surface spherical harmonics thus arise as particular cases of more general functions, 
viz., the ring-shaped functions. 

Let us go now to the radial equation. The introduction of (9) into (7) leads to 

^(^)-M ;+ ^)(' + ^ +i ) fl+2 ( E ^)^ ' (15) 

which is reminiscent of the radial equation for the hydrogen atom except that the orbital 
quantum number I is replaced here by I + (l/2)(5i + 5 2 ). The solution of (15) for the discrete 
spectrum is 

R(r) = i^(r; 6 U 5 2 ) = C nl {5 u 5 2 )(er) l+h ^e-^ 1 F 1 (-n + I + 1, 21 + 5 1 + 5 2 + 2; er), (16) 
where n G N — {0}. In Eq. (16), the normalization factor C n [(5i, 5 2 ) reads 



n (r r x 2(-0)* 1 

^nl{0l,0 2 ) 



( n+ fi±^) 2 r(2/ + 5 1 + 5 2 + 2)^ (n-l-1) 
and the parameter e is defined by 



T(n + I + 5 1 + 5 2 + 1) 



—^TT- (18) 



The eigenvalues E are then given by 



a 2 



E = E n = —-2, n = l,2,3,---, (19) 

2 (n + 

in agreement with Ref. (see also Refs. [9, 17]). Equations (16) to (19) can be specialized to 
the cases 5i = 5 2 = 5 (with 5 ^ or 5 = 0). In the limiting case 5 = 0, we recover the familiar 
results for hydrogenlike atoms. 

3 . Bi- Orthogonality of the Radial Wavefunctions 

The radial wavefunctions R n i [see Eq. (16)] satisfy the orthogonality relation 



Inn' 



/ Rn'i(r;6x,6 2 )Rni(r;6 1 ,6 2 )r 2 dr = 6 nn >, (20) 
Jo 



which is necessary for the normalization of the total wavefunction ^ n i m (r, 8, tp; 8%, 82) given by 
Eqs. (5), (8), and (16). In addition to the condition (20), we have also the following orthogonality 
condition 

f°° 2a 2 1 
Jw = R n v{r;8 1 ,8 2 )R n i{r;8 1 ,8 2 )dr = -3—. 8 W . (21) 

JO f n _|_ Si+5 2 y 21 + Oi + 6 2 + 1 

This new relation shall prove useful when dealing with the interbasis expansions in Section 5. 
The proof of Eq. (21) is as follows. 

In the integral appearing in (21), we replace R n i and R n i> by their expressions (16). Then, 
we take the hypergeometric function iF^—n + I + 1, 21 + 81 + 82 + 2; er) in (16) as the finite 
sum 

lFl (-n + l + l,2l + 8 1 + 82 + 2; er) = "j? J7] + !t ^ ^ 

(2/ + 81 + d 2 + 2) s s\ 



and we perform the integration term by term with the help of the formula [|3£ 

I> + 1 



By using 



we arrive at 



J™ e-^x\F x (a, r , kx) dx = ^± ] 2 F 1 I a,u + 1, 7; j J - (23) 



„ , , . T(c)T(c — a — b) , . 

2 F 1 (a 7 6, c; 1 = ' 1 -(, 24 

T(c — a)T(c — 0) 



2a 2 T(l + V + 8 1 + 6 2 + l) 



1 T(n + l + 8 1 + 8 2 + l) 



( n+ 5i±&) 3 r(2/ + 8 1 + 82 + 2) \(n-l-i)\(n-l'-i)\r(n + l' + 8 1 + 8 2 + l) 

"v^ 1 (~n + / + l) s (l + l' + 8 l + 8 2 + l) s T(n-l- s) 

s\ (2l + 8 1 + 8 2 + 2) s T(l'-l-s + lY [ ' 



By introducing the formula |39 



T(z-n) v 1 r(-z + l) 
into (25), the sum over s can be expressed in terms of the 2R1 Gauss hypergeometric function 
of argument 1. We thus obtain 

2a 2 1 



J w 



(n + 



S1M2.Y I + V + 81 + 82 + I 



n 



T{n + l + 5 1 + S 2 + l) 



\(n-l'-l)\ T(n + l' + 5 1 + 5 2 + l) T(l + l)T(l' -1 + 1)' 
Equation (21) then easily follows from (27) since [T(l - V + 1)T(V - I + l)]" 1 = 5 W . 



(27) 



The result provided by Eq. (21) generalizes the one for the hydrogen atom |39[]. Indeed, 
orthogonality properties similar to (21) hold for the hydrogen atom system and the harmonic 
oscillator system in /-dimensional spaces (/ > 2) |39[] . Such unusual orthogonality properties 
are connected to the accidental degeneracies of the energy spectra for these systems. In this 
connection, the property (21) is a consequence of the accidental degeneracy (with respect to the 
orbital quantum number /) of the discrete energy spectrum for the generalized Kepler-Coulomb 
system under consideration ; this point is further studied in the appendix. 



4 . Parabolic Basis 



In the parabolic coordinates 



x = ^JJL^>cos^p } y = y/JIUsm <p, z=—(ji — u), 0<fi,u<oo, 0<(f<2n, 



(28) 



the potential V reads 



2 2 2 
V = a h ci— r + c 2 - 



fM + U /i(/i + z/) v(fi + v) 
By looking for a solution of the Schrodinger equation (4) for this potential in the form 



(29) 



(30) 



we obtain the two coupled equations 

— (p— 1 + 
dfj, \ d\i) 



E (\m\+Si 



2^ 



dv \ dv I 



4/i 

E (|m|+5 2 ) 2 



-v — 



+ B 1 
+ B 2 



fi = 0, 



/ 2 = 0, 



(31) 



(32) 



2 Av 

where the two separation constants B\ and B 2 obey B\ + B 2 = —a. (In the whole paper, we 
use \& to denote the total wavefunction whatsoever the coordinate system is. The wavefuctions 
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\1/ in spherical, parabolic and prolate spheroidal coordinates are then distinguished by the 
corresponding quantum numbers.) By solving (31) and (32), we get the normalized wavefunction 



E 2 e imtp 



*(//, u, ip) = ^ ni „ 2m (/i, v, <p; 5 h 8 2 ) = r — f nim (p,; Sx)f n2m (u; £2) -7==, (33) 



where 



T(rii + \m\ +5i + l) 



T(\m\ + 5i + rii 



|m|+fj 

|t<J 2 e"t* 1 iFj (-Tit, H + Si + 1; , (34) 
with i — 1, 2 (ti = jj, and t 2 = v). In Eq. (34), we have 

n, = - H+ / < + 1 + ^ (* = 1 ' 2 )- ( 35 ) 

Here again m 6 Z in order to ensure the univaluedness of ^ nin2m . In addition, n, G N (for 
i = 1,2) in order that \E f niri2m be in L 2 (R 3 ). Therefore, the quantized values of the energy E 
are given by (19) where now the quantum number n is n = ri\ + n 2 + \m\ + 1, a number that 
parallels the principal quantum number of the hydrogen atom in parabolic coordinates. 

The results (33) to (35) agree with the corresponding ones of Ref. |19| derived by making 
use of the Kustaanheimo-Stiefel transformation (see also Ref. PI). 



Following Kibler and Campigotto [|n|, we can obtain a further integral of motion besides 
E. In parabolic coordinates, this integral corresponds to the hermitean operator 



x 2 



+ V 

The eigenvalues of X are 



d_( d_\ _ d_ ( cT 
dv I dv J V dfi \ dfi 



H-vd 2 2c x v 2c 2 /j, \i - v 

tt^7 "I 7 r 7 - — a . (do) 

2\iv dip 2 /i(/i + v) v[ji + v) fi + v 



ni -n 2 + ^ 



B 1 -B 2 = -a— r a 2 . (37) 

In Cartesian coordinates, the operator X can be rewriten as 



^ f d 2 ^ d 2 \ d 2 d 2 d z ^ r — z r + z 

\dx 2 dy 2 J dxdz dydz dz r r(r + z) r(r — z) J 



so that it immediately follows that X is connected to the ^-component A z of the Laplace- Runge- 
Lenz-Pauli vector via 

T — 7 V A- 7 

X = A z + Cl ^— - c 2 ^- (39) 
r(r + z) r(r — z) 

and coincides with A z when (3 = 7 = 0. 

5 . Interbasis Expansion between Parabolic and 

Spherical Bases 

The connection between the spherical (r, 9, ip) and parabolic (/i, u, if) coordinates is 

/j, — r(l — cos#), v = r(l + cos#), (^(parabolic) = y? (spherical). (40) 

For a given subspace corresponding to a fixed value of E n , we can write the parabolic wave- 
function (30) in terms of the spherical wavefunctions (5) as 

n-l 

is, ip; 5 U 62) = Yl W nm 2 m(^ S 2 )^nim(r, 9, f, 5 U 5 2 ) . (41) 

l=\m\ 



By virtue of Eq. (40), the left-hand side of (41) can be rewritten in spherical coordinates. Then, 
by substituting 9 = in the so-obtained equation and by taking into account that 



P M {1) = (^±l)n > 



(42) 

n\ 

we get an equation which depends only on the variable r. Thus, we can use the orthonormality 
relation (21) on the orbital quantum numbers I. This yields 

m-\m\ 

" 1 " jml1 ' 2J T{\m\+5 1 + l)T{2l + 5 1 + 5 2 + 2) ni " 2 nri1 ' 1 ; 

where 



E%r» = y/(2l + S 1 + S 2 + l)(l-\m\)\T(l + S 1 + l) 

1/2 





m 


+ Si + l)r(n 2 + 


m 


+ 5 2 + i)r 


fa 


+ / + 6i + 5 2 + 1) 


(m)!(n 2 )!(n-/-l)!r(/ + 5 ; 


> + i)r(/ + 


m 


+ 5 ± + 8 2 + 1) 



(44) 
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and 



K 



I'm 
nni 



X 



Z+|m|+5i+(52 



e~ x iF 1 (-n 1 , \m\ +S 1 + 1; x) iF^-n + l + l,2l + 5 1 + S 2 + 2;x) dx. (45) 



To calculate the integral , it is sufficient to write the confluent hypergeometric function 
n%, \m\ + 81 + 1; x) as a series, to integrate according to (23) and to use the formula (24) 
for the summation of the Gauss hypergeometric function 2 F 1 (a, b, c; 1). We thus obtain 



K 



1 in 

nni 



[n — 


m 


-l)!r(2Z + 5 1 + 5 2 + 2)r(Z + 


m 


+ 81 + 8 2 + 1) 


(I- 


m\)\T(n + l + 5i + 5 2 + l) 



p \ —ni, —l + \m\,l+ \m\ + 81 + 8 2 + 1 
' 3 2 \ \m\ + 81 + 1, -n + \m\ + 1 

The introduction of (46) into (43) gives 



W, 



(61, 8 2 



m — m 



(n — 


m 




-i)!r(z + 


m 


+ 5i + 8 2 + 1) 


(i- 


m|)!r(| 


m 


+ St + l)r(n + / + 5i + 5 2 + 1) 



7 

-1)- 



-ni, —I + ^ + \ m \ + 81 + 8 2 + 1 



3 2 y \m\ +81 + I, —n + \m\ + 1 
and owing to (44) we end up with 



l\E. 



Irn 
nin 2 



\ 



(21 + 8i + 8 2 + i)r(m + 


m 


+ 5i + i)r(n 2 + 


m 


+ 8 2 + 1) 


(ni)!(n 2 )!(n - 


/-!)!(/- 


m )! 



(n — 


m 


-1)! 


r{i + 8i + i)v{i + 


m 


+ 81 + 8 2 + 1) 


r( m 


+ 5i + l)\ 


T{1 + 8 2 + l)r(n + Z + 5i + 8 2 + 1) 



m - 

•(-1)- 

p [ -ni, —1+ \m\,l+ \m\ + 5i + 8 2 + 1 
' 3 2 \ \m\ + 81 + 1, —n + \m\ + 1 

that constitutes a closed form expression for the interbasis coefficients. 



(46) 



(47) 



(48) 



The next step is to show that the interbasis coefficients (48) are indeed a continuation on 
the real line of the Clebsch-Gordan coefficients for the group SU(2). It is known that the 
Clebsch-Gordan coefficient C^. 6 * = (aabj3\abcy) can be written as |37j 



aa;b(3 



{-l) a - a 8^ a+p {a + b - 7 )!(6 + c - a)\ 
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(2c + l)(a + a)!(c + 7 )! 



(a - a)\(b - P)\(b + (3)l(c - j)\(a + b + c+ l)!(a + b - c)!(a - b + c)!(6 - a + c)\ 



1/2 



3-r 2 



-a — 6 — c — 1, —a + a, — c + 7 
-a — 6 + 7, — b — c + a 



In view of the formula [HO 



3^ 2 <| ./ 



t'A-N-t 



(t + s) N j s ,t'-s',-N 

3^ 2 S -/ 



1 , 



equation (49) can be rewritten in the form 

(2c + 1)(6 - a + c)!(a + a)!(6 + /3)!(c + 7)! 



^aa;bl3 



-D a - a s, 



a - a)!(6 - /?)!(c - j)\(a + b- c)!(a - b + c)!(a + 6 + c + 1)! 

il 



(a + b- 7)! f -a + a,c + 7 + l,-c + 7 
(6 -a + 7)! 3 2 \7 -a-&,&-a + 7 + l 



By comparing Eqs. (51) and (48), we finally obtain 



Wi 



, m-\m\ 1 I S l + 6 2 | m | + 21 

HH 5^ 2 >l'"l^ 2 

^- / n+(52 — 1 I +712 — n l "t"^2 1 |m|+ni — 7i2+<5^ 



(49) 



(50) 



1/2 



(51) 



(52) 



Equation (52) proves that the coefficients for the expansion of the parabolic basis in terms of 
the spherical basis are nothing but the analytic continuation, for real values of their arguments, 
of the SU(2) Clebsch-Gordan coefficients. 
The inverse of Eq. (41), namely 



n— \m\ — 1 



(ji, u,(p;5 h 5 2 



(53) 



is an immediate consequence of the orthonormality property of the SU(2) Clebsch-Gordan 
coefficients. The expansion coefficients in (53) are thus given by 



m^mj 7 + £l+j2 | | m | + £l+j2 

" + ' 5 2 2 ~ 1 , " + a 2 " 1 -ni; " + 2 1 ,m+M- g=% -1 



(54) 



and may be expressed in terms of the 3_F 2 function through (49) or (51). 

To close this section, it should be mentioned that (52) and (54) generalize the well-known 
result, corresponding to 5i — 82 — 0, for the interbasis expansion between parabolic and 
spherical bases obtained in Refs. [41-44] in the case of the hydrogen atom. Furthermore, by 
taking 8\ = 5 2 7^ in (52) and (54), we recover our former result [II] for the Hartmann system. 
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6 . Prolate Spheroidal Basis 



6.1. Separation in Prolate Spheroidal Coordinates 
We now pass to the prolate spheroidal coordinates 



R 



x 



R 



R 



o v (e 2 -l)(l-^)cos^ 2 / = - v /(e 2 -l)(l-r ? 2 )sin^, z = -(^ + l), (55) 



1 < £ < 00, 



T < 77 < 1, 



< < 2tt, 



where i? is the interfocus distance. In the limits where R — > and i? — > oo, the prolate 
spheroidal coordinates give back the spherical coordinates and the parabolic coordinates, re- 
spectively |27|, In the system of prolate spheroidal coordinates, the potential V can be 
written as 



V = a 



+ 



l 



cr 



+ c 2 - 



1 



(56) 



The Schrodinger equation (4) for the potential (56) is separable in prolate spheroidal coordi- 
nates. As a point of fact, by looking for a solution of this equation in the form 



(57) 



we obtain the two ordinary differential equations 



d r , 2 n d , (H+^) 2 (H+£ 2 ) 2 0nm E# 2 



-\(R)i/)i, (58) 



dr/ 77 > dr] 2(1 + 7]) 2(1-7]) 



2aR7] + 



ER 2 



[i-v 2 



^ = -A(i?)V 2 , (59) 



where X(R) is a separation constant in prolate spheroidal coordinates. By eliminating the 
energy E from Eqs. (58) and (59), we produce the operator 



A 



£ 2 -7f 



(Z 2 - 1)(1 - 7] 2 ) dtp 2 



_ aR ^L + 1 + 2 AZ + v) 2 + (Z-m-v) + 2c2 (e + V) 2 ~ it + 1)(1 + r?) 



(e + ^)(e + l)(l + 77) 



(60) 
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the eigenvalues of which are X(R) and the eigenfunctions of which are ipiiQfaiv) ■ The signif- 
icance of the (self-adjoint) operator A can be found by switching to Cartesian coordinates. A 
long calculation gives 

A = M + RX, (61) 
where X is the constant of motion (38) and M is the following constant of motion 

M = L 2 + Cl ^^ + C2 -1^, (62) 
cos z 2 sin 2 

the operator L 2 being the square of the angular momentum operator. The prolate spheroidal 
wavefunctions ipi and ip 2 could be obtained by solving (58) and (59). However, it is more 
economical to proceed in the following way that shall give us, at the same time, the global 
wavefunction (57), i.e., rj, ip; R; 5±, 8 2 ), and the interbasis expansion coefficients. 

6.2. Interbasis Expansions for the Prolate Spheroidal Wavefunctions 

From what preceds, we have three sets of commuting operators, viz., {H, L z , M}, {H, L z , X}, 
and {H, L Z ,A} corresponding to the spherical, parabolic, and prolate spheroidal coordinates, 
respectively. (The operators L z and H are the z-component of the angular momentum and the 
Hamiltonian for the generalized Kepler- Coulomb system, respectively.) In particular, we have 

M* nlm (r, 6, p; 5 U 5 2 ) = (l + (l + + X ) 9 > ^ 5 ^ ^ ' ^ 

X^ nin2m (^ v, p; 5i, 5 2 ) = -a ■ h+h 2 *„ in2 m(^, v, p; 5 U S 2 ), (64) 

n+ 2 

and 

A^ nqm (£, 7], ip; R; 8 ± , 5 2 ) = X q (R)^ nqm (^, V, V 5 ; R\ <*i, fa) (65) 

for the spherical, parabolic, and prolate spheroidal bases, respectively. In Eq. (65), the index q 
labels the eigenvalues of the operator A and varies in the range < q < n — \m\ — 1. We are 
now ready to deal with the interbasis expansions 

n-1 

^nqm(^ V, <P\ R\ fa, fa) = J2 T nqm( R 'i fa, fa)^nlm(r, 0, p; 5 U 5 2 ), (66) 

l=\m\ 
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ft— m — 1 



y nqm {£, V, W R; Si, 6 2 ) = U nqm( R '> S l, ^)^n 1 n 2 m(/i, V, <f] 8 X , S 2 ) 



(67) 



ni=0 



for the prolate spheroidal basis in terms of the spherical and parabolic bases. [Equation (66) 
was first considered by Coulson and Joseph [p8[] in the particular case 61 — 5 2 — 0.] 

First, we consider Eq. (66). Let the operator A act on both sides of (66). Then, by using 
Eqs. (61), (63), and (65) as well as the orthonormality property of the spherical basis, we find 
that 



X q (R) - I 



5i + 5 2 



1 + 



Sx + S 2 



+ 1 



n-l 



Tn qm {R] Si, 5 2 ) - R 53 T l nqm (R;8i,5 2 )(X) w , 



l'=\m\ 



(68) 



where 



(XV = Kim(r, 6, ip; 5 U 5 2 )Xm nVm {r, 6, <p; 5 U 5 2 ) r 2 sin 9drd9d<p. (69) 

Jo Jo Jo 

The calculation of the matrix element (X)hi can be done by expanding the spherical wave- 
functions in (69) in terms of parabolic wavefunctions [see Eq. (53)] and by making use of the 
eigenvalue equation for X [see Eq. (64)]. This leads to 



a 



n— ra — 1 



n 



£ (2n 1 -n+\m\ + 5 -^ + l)w:USi,S2)W: i ) m (Si,S2). (70) 

2 ni=0 V Z / 



Then, by using Eq. (54) together with the recursion relation j37| 

4c 2 (2c+ l)(2c- 1) 



nl/2 



(c + 7)(c - 7)(-o + b + c)(a - b + c)(a + b - c + l)(o + & + c + 1) 



(c - 7 - l)(c + 7 - l)(-a + & + c- l)(a - 6 + c- l)(a + b - c + 2) (a + b + c) 

4(c- l) 2 (2c-3)(2c- 1) 

(a - /3)c(c - 1) - 7 a(a + 1) + 76(6 + 1) 



-,1/2 



c; 



2cfc- r 



c-2,7 

aa;b/3 



(71) 



and the orthonormality condition 



E 

«,/3 



°aa;6/3 u oc«;6/3 — c / c 7 / 7 , 



(72) 
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we find that {X)ip is given by 
2a 



(X)li' 



2n + <Ji + 8 2 



a(2\m\ + 8 1 + 5 2 )(5 1 - 8 2 ) 
(2l + 5 1 + 5 2 )(2l + 5 1 + 5 2 + 2) 



5m, (73) 



where 



A 



21 + 5! + 6 2 



(/ - \m\)(l + \m\ +St + 8 2 ){l + Si)(l + 5 2 )(n - l)(n + I + 5 X + 8 2 ) 
(2l + 5 1 + 5 2 -l)(2l + 5 1 + 8 2 + l) 



1/2 



(74) 



Now by introducing (73) into (68), we get the following three-term recursion relation for the 



coefficient T l m 



\+ 6 -^] (l+ 6 -^ + l 



aR(2\m\+ 5 1 + 5 2 )(5 1 - 5 2 ) 



+ 



2 
aR 



, (2l + 5 1 + 5 2 )(2l + 5 1 + 5 2 + 2) " q 
A l nmT^(R; 81,82) + A l nm T l nq ^(R] 8 X , 8 2 



X q (R) 



T nqm( R 'i ^1, 8 2 ) 



0. 



(75) 



2n + Sx + 5 2 

The recursion relation (75) provides us with a system of n — \m\ linear homogeneous equations 
which can be solved by taking into account the normalization condition 



n-l 



E \T'(R;5 1 ,52)\ 2 = l. 



(76) 



l=\m\ 



The eigenvalues X q (R) of the operator A then follow from the vanishing of the determinant for 
the latter system. 

Second, let us concentrate on the expansion (67) of the prolate spheroidal basis in terms of 
the parabolic basis. By employing a technique similar to the one used for deriving Eq. (68), we 
get 



X q (R) + aR 



ni -n 2 + 



n + 



5i+8 2 



where 



n— |m| — 1 

U% m (R;6 1 ,6 2 ) = £ U< m (R;8 1 ,8 2 )(M) nin{ , (77) 

n' 1= 



rcc roo rZir u J- y 

(M) nin /=jf J J ^ in2m (/i,z/,^;5 1 ,5 2 )M^ n; ^ m ( y u,i/,^;5 1 ,5 2 )^^rf/idz/^. (78) 

The matrix elements (M) nin ' can be calculated in the same way as (X)u> except that now we 
must use the relation Wi 



l(b-a + c)(a-b + c+ l)] 1/2 CZ., b , = [{a-a + l){b- /?)] 1/2 ^ : 



i+l/2,a-l/2;6-l/2,/3+l/2 
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+ [(a + a + 1)(6 + /9)] 1/2 ^ lAQ+1/2;fe _ 1/2j/3 _ 1/2 
and the ort honor mality condition 

E C al;b/3 C ll>;bf3> = <W<W, 
c,7 

instead of Eqs. (71) and (72). This yields the matrix element 



(79) 



(80) 



(n — n-i — \m\ - l)(n x + 1) + (n x + |m| + S 1 )(n — n 1 + S 2 ) + -5 2 -2) 



-[(n - ni - \m\ - l)(m + l)(ni + |m| + 5i + - ni + o~ 2 - 1)] 1/2 o^.m+i 

- [m(n - ni - |m|)(ni + |m| + o\)(n - m + S 2 )] 1/2 oV^m-i- (8 1 ) 
Finally, the introduction of (81) into (77) leads to the three-term recursion relation 



(n — ni — \m\ - l)(ni + l)(m + \m\ + 5i)(n — ni + 8 2 ) + -(o~i - S 2 )(5 1 - 5 2 - 2) 



—aR 



2m - n+ \m\ + 1 + 
n + 



<5i+<5 2 
2 



- A,(i2) 



c/- m (^;5i,5 2 ) 



[(n -m- \m\ - l)(m + l)(m + |m| + <$! + l)(n - m + 5 2 - 1)] 1/2 [/^(i?; S u 5 2 ) 

+ [ ni ( n _ ni _ | m |)( ni + | m | + Sl )(n - m + 5 2 )] 1/2 ^ 1 m 1 ( J R; 5i, <f 2 ) (82) 

for the expansion coefficients U^ m (R;5i,5 2 ). This relation can be iterated by taking account 
of the normalization condition 



n— \m\— 1 

E |t^m(i2;<yi,«fe)i 2 = i. 

m=0 



(83) 



Here again, the eigenvalues \ q (R) may be obtained by solving a system of n — |m| linear 
homogeneous equations. 

In the case 5i — 8 2 — 0, from Eqs. (75) and (82) we obtain three-term recursion relations for 
the coefficients of interbasis expansions of the prolate spheroidal basis in spherical and parabolic 
bases for the hydrogen atom ; these coefficients were calculated in Refs. [28, 44, 45]. 
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Finally, it should be noted that the following four limits 

Jim U% m (R; 5 U 5 2 ) = 5 niq , Km U% m (R; 5 U 5 2 ) = W% m (6 u 5 2 ), (84) 

|m T l nqm {R- 5 U 5 2 ) = 6 k , Jim T l nqm (R; 5 U 5 2 ) = W^S,, 5 2 ) (85) 
furnish a useful means for checking the calculations presented in Sections 5 and 6. 

7 . Concluding Remarks 

One of the main results of this paper concerns the construction of the prolate spheroidal 
basis for the generalized Kepler-Coulomb system (corresponding to 5\ ^ and 5 2 ^ 0) as a 
superposition of either the spherical basis or the parabolic basis. In this respect, our work is a 
continuation of the pioneer work by Coulson and Joseph p8| on the ordinary Kepler-Coulomb 
system (corresponding to 8\ — and 8 2 — 0). We have obtained three-term recursion relations 
for the coefficients of the interbasis expansions of the prolate spheroidal basis in terms of the 
spherical basis, on one hand, and of the parabolic basis, on the other hand. These recursion 
relations can be easily implemented by using symbolic programming languages (like MAPLE, 
MATHEMATICA, etc.) especially, in the framework of perturbation theory, for small or large 
values of the spheroidal separation constant R. 

We have also calculated the coefficients of interbasis expansions between the spherical basis 
and the parabolic basis. As we could expect from the limiting case 5\ = 5 2 = 0, the latter 
coefficients are Clebsch-Gordan coefficients for the special unitary group SU(2) modulo an 
analytic continuation to real values of their arguments. 

It was realized long time ago that prolate spheroidal coordinates may be useful in quantum 
chemistry even in situations exhibiting spherical symmetry [28, 29]. This work shows the rele- 
vance of prolate spheroidal coordinates in axial symmetry. This is indeed the kind of symmetry 
that is encountered in the theory of the chemical bond. 

We close this paper with a remark of interest for mathematical physics (and mathematical 
chemistry). Most of the 5±- and ^-dependent expressions (e.g., the ring-shaped functions) in 

18 



this work present nice properties with respect to the arguments 5 1 and 5 2 - It would be worth 
to study some of these expressions in connection with the theory of orthogonal polynomials. 
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Appendix: An Alternative Proof for the 
Bi-Orthogonality of the Radial Wavefunctions 

We start from two copies of Eq. (15): one for the pair (E = E n , R = R n i) and the other one 
for (E = E n /,R = Rn'i'). After multiplication of the first copy (respectively second copy) by 
(respectively i? n ;) and integration, with the measure r 2 <ir, on the half real line, we obtain 

(l'-l)(l + l' + S 1 + S 2 + l) / R n n>(r;6 1 ,6 2 )R nl (r;6 1 ,6 2 )dr 

Jo 

= 2{E n -E n .) R n , l ,(r;6 1 ,6 2 )R nl {r;6 1 ,5 2 )r 2 dr. (86) 
Jo 

From the latter equation, it readily follows that 

PCX) 

/ R n > l (r;6 1 ,6 2 )R nl (r;6 1 ,6 2 ydr = for n' ^ n (87) 
Jo 

and 

/•oo 

/ Rra>{r;6 1 ,6 2 )R nl (r;6 1 ,6 2 )dr = iorl'^l. (88) 
This completes the proof of the two orthogonality relations 

Inn' ~ vnn' 

and J w ~ <5;/'. It should 

be emphasized that Eq. (87) turns out to be a direct consequence of the independence of the 
energy E on the quantum number I. 
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